00 

o 

O 
(N 



Massive particle production to NNLO in QCD 
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We discuss the recent derivation of the one-loop squared virtual QCD corrections to the W boson pair production 
in the quark-anti-quark-annihilation channel in the limit where all kinematical invariants are large compared to 
the mass of the W boson. In particular, we elaborate on the combined use of the helicity matrix formalism with 
the Mellin-Barnes representations technique. 



1. Introduction 

The Large Hadron Collider (LHC) is expected 
to have a huge impact on particle physics phe- 
nomenology. Most of the processes which will be 
studied at the LHC need to be calculated at least 
to next-to-leading order (NLO) in QCD whereas 
there are some for which a theoretical prediction 
is needed to next-next-to-leading order (NNLO). 
Electroweak gauge boson pair production falls 
into the latter category. One of the reasons is 
that the increase of the centre-of-mass energy at 
the LHC with respect to the Tevatron from 1.96 
TeV to 14 TeV will result in a huge boost of the 
available data. 

The importance of hadronic W-pair produc- 
tion is two-fold. Firstly, it is a process which 
allows the measurement of the vector boson tri- 
linear couplings and therefore a comparison with 
the Standard Model (SM) predictions. Most 
attempts to model New Physics, such as Su- 
persymmetry and Extra-dimensions in all varia- 
tions, should be able to explain any deviations by 
consistently adjusting the anomalous couplings 
and/or by incorporating decays of new particles 
into vector boson pairs [1,2]. 

Secondly, hadronic W pair production is im- 
portant for investigations of the nature of the 
Electroweak symmetry mechanism by contribut- 
ing the dominant background for the Higgs boson 
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mediated process (see Refs. [3-15]), 

pp ^ H ^ W*W* ^ IvT'v' , 

in the Higgs mass range between 140 GeV < 
Mh < 180 GeV [16]. 

The interest in hadronic W pair production is 
well displayed by the fact that the Born cross sec- 
tion was calculated some thirty years ago [17]. 
The NLO QCD corrections were computed in the 
90's and seen to contribute a 30% [18-22]. Next, 
soft gluon resummation effects were considered 
in Ref. [23] whereas massless fermion-boson scat- 
tering was studied at NNLO in Ref. [24]. The 
first steps towards a complete NNLO study were 
taken with the computation of the NNLO two- 
loop [25-27], as well as the one- loop squared [28] 
virtual corrections in a high energy expansion, 
< s, t, u. 

The methods used in [28] , differ somehow from 
the ones employed in [26] though, both are a con- 
tinuation of the techniques used before in [29-33] , 
The difference lies mainly in the fact that for 
the one-loop squared corrections the helicity ma- 
trix formalism was additionally used to reduce 
the problem to a small set of integrals, which in 
turn were treated with Mellin-Barnes (MB) rep- 
resentations [34,35]. The latter were constructed 
by means of the MBrepresentation package [36] 
and then analytically continued in the number of 
space-time dimensions _D = 4 — 2e using the MB 
package [37]. After the asymptotic expansion in 
the mass parameter, contours were closed and in- 
tegrals finally rcsummcd cither with the help of 
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XSummer [38] or the PSLQ algorithm [39]. 

Here, we are going to give more details on how 
the two methods of the helicity formalism and MB 
representations were combined for the derivation 
of the result in [28]. We provide, as an exam- 
ple, the coefBcient of a certain helicity matrix el- 
ement for the one-loop amplitude in the high en- 
ergy limit. This result is in closed analytic form 
expressed through harmonic polylogarithms and 
transcendental constants. 

2. The Calculation 

We shall introduce here part of the notation 
used in [28]. The charged vector-boson produc- 
tion in the leading partonic scattering process cor- 
responds to 

qj{pi) + qj{p2) -^ W^{p3,m) + W+{pi,m),{l) 

where pi denote the quark and W momenta, m is 
the mass of the W boson and j is a flavour index. 
Here we are considering down-type quark scatter- 
ing. Energy-momentum conservation implies 
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We consider the scattering amplitude A4 for the 
process (H]) at fixed values of the external parton 



momenta pi, thus pf 



pI 



and pI — pj 



m . The amplitude A4 may be written as a series 
expansion in the strong coupling Og, 



\M) 



|A^(o)) 



27r 



|A^(i)) 



2,71 



(3) 



For convenience, we define the function 
A(t,m,s,t,^) for the squared amplitudes 
summed over spins and colors as 



spin, color 



w+ + w-)\'^ = 

A{e,m,s,t,^) . 



(4) 



^ is a function of the Mandelstam variables s, t 
and u given by s = {pi +P2)^, t = (pi -pzY ~m? 
and u = {pi — pi)'^ — rn? and has a perturbative 



expansion similar to Eq. ([3]) 
A{e,m,s,t,iJL) = 



^(0)+[^]^(l) 

27r 



27r 



^(2) + 0{al) 



(5) 



In terms of the amplitudes the expansion coeffi- 
cients in Eq. (O may be expressed as 



+ (Al(2)|Al(o)) 



(6) 

(7) 

(8) 



As already mentioned, in order to compute 
(A^(i)|A^(i)) we used the helicity matrix formal- 
ism, namely we expressed the result in terms of 
hehcity amplitudes, A^^(Ai, A2, s,i). The quark 
and the anti-quark have opposite helicities in the 
centre-of-mass system so one helicity label above, 
g — ±1, suffices. Ai and A2 stand for the hehcities 
of the W"*" and W~ respectively. 

Starting from the one-loop amplitude, the ini- 
tial expression can be rearranged as 
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Y^Cf{s,t,u,m)M 



g 

i ' 



(9) 
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where the Ci are coefficients and A^f are helicity 
matrix elements and g = ±. The ten helicity 
matrix elements have been taken as defined in 
Ref. [40] (see also [41]): 
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Ml 
Ml 
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Ml 



viP2)(/iil^ -p/2)el2Vgu{pi), 

v(p2)lj'3'PgU{pi)ei -62, 
V{p2) €/i_Vg U{pi) 62 -ps, 

-v{p2) f^Vg u{pi) ei • P4 , 

V{p2) e/lVg U{pi) €2 -pi, 

-v{p2) e/iVg u{pi) ei ■p2, 
HP2) T^zPg u{pi) ei • _P2 £2 ■ Pi 
v{p2) J^zPg u{pi) ei • _P2 £2 • P3 
v{p2)l^zPgu[pi)ti •P4e2 -Pi 
v(p2)lf3'Pg u{pi) ei • _P4 £2 • Pa 



(10) 
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where Vg = V± 



i^^. All colour indices as 



well as the arguments of the polarization vectors, 
ei(p3, Ai) and £2(^41 ^2), have been suppressed. 

Let us have a look at the CQ{s,t,u,m) coeffi- 
cient of the M.Q matrix element for the one-loop 
amplitude. Typically, C^ {s,t,u^m) is given by 

Cq {s^t^u, m) = 2^ Co,i(s, t, u, m)Ii{s, t, u, m; fi^) (11) 

i 

where cg ,i(s, t, u, m) are polynomials in the kine- 
matical variables and Ii{s,t,u,m] fi^) are one- 
loop integrals. For example, one of the one-loop 
integrals appearing in Eq. [TT] is 

k.pi 



I 



d'^k 



k'^{k+plY{k+pAy 



(12) 



After feeding this into the MBrepresentation 
package, one gets, from the reduction of the orig- 
inal tensor structure into scalar objects, the fol- 
lowing two terms: 
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xr(-e)r(l-e-zi)r(-zi) 




xr(i + zi)r(n-e-t-.2i) 



(13) 



x(r(2-2e))- 



(14) 



where of course / — if^^ + /^^ . One then needs 
to perform an asymptotic expansion in the mass 
parameter and finally resum the MB integrals. 
This allows to compute the coefficients Cf in 
a closed analytic form in the high energy limit. 
As an example, we present here the result up to 
0{e^) for Cq [Cp the color factor and g^^ the 
quark-W coupling) 



Cq — Cp gy^i^ 



+ - 



^' \ 
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-j-ZTT 
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1-x 



T3 



3(1 -x) 



10 3,2 

X 1-x I y 
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1 f 7C3 



8(1 -x) 3 Vl-a; 



27 



81,2(2;) 



1 — x J 1 — X 

L^ _ Li2(x) _ 5hy 
2(1 — x) 1 — X X 

Stt^ U2{x)tt'^ 



4(1 -x) 
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737r4 L^ 
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1 /20 3 
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57r2 3 
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1 f 7C3 36_^| ^ fl 1^20 
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2 V 1 - a; 1-x 
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+i'K e 



H 



1 — X 1 — X 

5Ll 



6(1 -a;) 2x 



4(1 -x) X J " 
7C3 , 5Li2(x) 



Li2(a;)Lj^ 
1-a; 



3(1 — a;) X 

Li3(a;) 81,2(3;) 



1 — 3; 1 — a; 
where we have defined 



(15) 



X = 
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— , y = - 
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-- , ms = 

s s 


(16) 


and 








-Lm = 


-- Log(ms) , 


Ly — Log (1 — 3;) . 


(17) 
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Similarly, one can compute the Cf coefficients 
for all the helicity matrix elements Adf. The com- 
plex parts of the coefficients are given explicitly 
as can be seen in Eq. [T5] which means that ob- 
taining the complex conjugate expressions for the 
tree-level, |A^*-°^), and the one loop, \^A^^^'), am- 
plitudes is trivial. The following step would be 
the contraction (A^*^^)|A^(^)). As an easy test, 
we have checked (A1*^°^|A1*^^^) derived with this 
method against the results provided in [18,19]. 

In order to check (A^*^^)|A^^^)) we have used a 
more involved test of the infrared structure ac- 
cording to the Catani prediction [42]. In the 
case of one-loop QCD amplitudes, their poles in 
e can be expressed as a universal combination of 
the tree amplitude and a colour-charge operator 
P^^(e). The generic form of T"'^'(e) was found by 
Catani and Seymour [43]. 

3. Conclusions 

We have discussed some details of the computa- 
tion of the one- loop squared NNLO QCD virtual 
corrections for the process qq — > W~^ W~ in the 
limit of small vector boson mass. Our main result 
was presented in [28] . This was a second step to- 
wards the complete evaluation of the virtual cor- 
rections. In a forthcoming publication, we will de- 
rive a series expansion in the mass and integrate 
the result numerically to recover the full mass de- 
pendence, similarly to what has been done in [44] . 

To complete the NNLO project one still needs 
to consider 2^3 real-virtual contributions and 
2^4 real ones. The real-virtual corrections are 
known from the NLO studies on WW + jet pro- 
duction in [45,46]. 
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